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Abstract 



We develop a method for computing form factors of local operators in the framework 
of Sklyanin's separation of variables (SOV) approach to quantum integrable systems. 
For that purpose, we consider the sine-Gordon model on a finite lattice and in finite 
dimensional cyclic representations as our main example. We first build our two central 
tools for computing matrix elements of local operators, namely, a generic determinant 
formula for the scalar products of states in the SOV framework and the reconstruction 
of local fields in terms of the separate variables. The general form factors are then 
obtained as sums of determinants of finite dimensional matrices, their matrix elements 
being given as weighted sums running over the separate variables and involving the 
Baxter Q-operator eigenvalues. 
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1 Introduction 



The computation of general matrix elements and correlation functions of local operators is one of the fun- 
damental problems of quantum field theory and statistical mechanics. These objects contain indeed the key 
dynamical and measurable quantities of the corresponding physical systems, see e.g. . In the integrable 
(low dimensional) situation 1151 4 1 in . thanks to the existence of powerful algebraic structures related to the 
Yang-Baxter algebra (see e.g. Ill2l - l20ll and references therein), significant progress towards their exact de- 
termination has been obtained in the last thirty years. Such results concern in particular models solvable by 



means of the algebraic Bethe ansatz like the XX Z Heisenberg spin chain ||5|-|9|]. They were first obtained at 
the free fermion point, namely for the Ising model and the Heisenberg chain (for anisotropy A equal to zero) 
112 1142411 . Going beyond such a free fermion case involves a deep use of the Yang-Baxter algebra. After his- 
torical attempts for the finite chain in the framework of the Bethe ansatz (see [18] and references therein) but 
leading in fact to implicit representations in terms of dual fields, explicit representations for form factors and 
correlation functions were first obtained directly in the infinite chain (in the massive regime) ll^Eill . The 
underlying quantum algebra structure was instrumental there together with a few assumptions on the repre- 
sentation of the Hamiltonian and of the local spin operators within the representation theory of this quantum 
affine algebra. Similar results for the disordered regime were then derived using the assumption that cor- 
relation functions should satisfy (/-deformed KZ equations in close analogy with the massive regime [26]. 
The derivation of these results (both for massive and massless regimes and in the presence of a magnetic 
field) was later obtained in the framework of the algebraic Bethe ansatz, starting from finite size systems, 
thanks to the resolution of the quantum inverse scattering problem B27I429I1 . Further investigations led to 
the extension of these results to the non-zero temperature case H30L bill and also to the case of non-trivial 
(integrable) boundary conditions 11321 13311 . The computation of physical (two-point) correlation functions 
required sophisticated summation techniques of the elementary blocks of correlation functions B34|-|37[| that 
ultimately led to the exact computation of their asymptotic behavior ll38l44lll. It is also worth mentioning 
that controlled numerical summation techniques combined with these exact results on form factors led to the 
determination of dynamical structure factors in very good agreement with actual neutron scattering experi- 
ments on magnetic crystals, see e.g. Il42|445ll . Another important line of research revealed powerful hidden 
Grassmann structures that could also be used in the analysis of the conformal limit of these models [|46j-|52j]. 
Finally, it was recently shown how to obtain the asymptotic behavior of the correlation functions for critical 
systems through their expansion in terms of form factors in the finite volume 115 3145 511 . 

For more sophisticated systems, like lattice integrable models related to higher rank algebras or for in- 
tegrable relativistic quantum field theories, the present state of the art is slightly less satisfactory despite 
considerable efforts. 

On the one hand, the bootstrap approach has provided great insights into the structure of the exact scattering 
matrices [56-62] and form factors of such integrable massive quantum field theories Il63l467ll. Perturbed 
conformal field theorjQ (CFT) 1 73 - 84 ] has also been used in this context together with a new understand- 
ing of the integrable structure of CFT I185l487ll . Several important attempts have also been made to use 
Sklyanin's separation of variable method (SOV) [88-94] in these more complicated algebraic situatio ns an d 
in particular in the case of infinite dimensional representations associated to the quantum fields [95-102]. 



See 168-72] and references therein for some literature on conformal field theories. 
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One should also mention that the new fermionic structures mentioned above in the case of the XX Z lattice 
model have been used recently to investigate the structure of matrix elements of the sine-Gordon model in 
the infinite volume limit B5 lL 15211 . 

On the other hand, although these advances lead to invaluable informations on these theories, a direct com- 
putation of the form factors and correlation functions starting from the description of their local fields and 
using the given Hamiltonian dynamics is still missing. The main reason is that although these approaches 
succeeded to describe more and more efficiently both the space of states of such models and the space of 
their form factors, the connection between the set of operators used to reach this goal and the local operators 
of the theory one is interested in has not been found yet. Consequently, the matrix elements and correlation 
functions of these local fields can at best be identified through indirect arguments. In other words, contrary 
to the above mentioned finite dimensional cases, the solution of the quantum inverse scattering problem is 
in general not known in the field theory situation; moreover, it appears very often in these more compli- 
cated models (either on the lattice or in the continuum) that the usual algebraic Bethe ansatz is no longer 
applicable due to the lack of a proper reference state, and that the SOV framework 1I88U90I1 has to be con- 
sidered instead. Although this beautiful method is quite general and powerful to describe the spectrum of 
these models, the problem of reconstructing the local operators of the corresponding theories in terms of the 
separate variables is in general still open (see however [99]). 

To explain in more detail the main issues of this problem, let us consider in particular the sine-Gordon model 
which will be the main example analyzed in this article. As we just recall, integrable massive field theories 
in infinite volume can be solved on-shel^ through the determination of their exact S -matrices [56-62] which 
completely characterize the particle dynamics. In this particle formulation of the theory a direct access to 
the local fields is missing and any information about them needs to be extracted from the particle dynamics. 
The monodromy properties and the singularity structure provide a set of functional equations [63-67] for the 
form factors of the local fields on asymptotic particle states. These form factor equations are uniquely fixed 
by the knowledge of the exact S -matrix and the space of their solutions is expected to coincide with that of 
the local opera tors of the theory; many results are known on the form factors of local fields, see for example 
1 78148 lh . il03l4ll3ll and references therein for some literature related to sine-Gordon model. However, it is 
worth pointing out that the form factor equations only contain information on symmetry data of the fields 
(like charges and spin) and after fixing them to some values, we are still left with an infinite dimensional 
space of form factor solutions which should correspond to the infinite dimensional space of local fields 
sharing these data. There is a large literature dedicated to the longstanding problem of the identification 
of the local fields in the scattering formulation of quantum field theories. Different methods have been 
introduced to address it; one important line of research is related to the description of massive integrable 
quantum field theories as (super-renormalizable) perturbations of conformal field theories by relevant local 
fields [73-84]. This characterization has led to the expectation that the perturbations do not change the 
structure of the local fields in this way leading to the attempt to classify the local field content of massive 
theories by that of the corresponding ultraviolet conformal field theories. The latter issue was first addressed 
in [80] in the simplest massive free theory, namely the Ising model. There a conjecture was introduced 
defining a correspondence between mild asymptotic behavior at high energy of form factors and chiral local 
fields. Such a conjecture was justified showing the isomorphism of the space of chiral local fields in the 



2 See for a review [83] and references therein. 
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massive and conformal models. The extension of the chiral isomorphism to interacting massive integrable 



theories was done in 111411 for s everal massive deformations of minimal conformal models, in 111 1511 for the 



sine-Gordon model and in 1 116 1 for all its reductions to unitary minimal model jE The prob lem to extend 



the classification also to non-chiral local fields was analyzed in a series of works in particular, 



for the massive Lee- Yang model, the first proof that the operator space determined by the particle dynamics 



coincides with that prescribed by conformal symmetry at criticality was given in II127I1 . While these are 
indubitably important results on the global structure of the operator space in massive theories it is worth 
pointing out that they do not lead to the full identification of particular local fieldsB In J93I1 a criterion has 
been introduced based on the quasi-classical characterization of the local fields; it has been fully described in 
the special cases of the restricted sine-Gordon model at the reflectionless points for chiral fields and verified 



on the basis of counting arguments [94]. 

This makes clear that, in the S-matrix formulation of massive quantum integrable theories in infinite vol- 
ume, the main open problem remains the absence of a direct reconstruction of the local fields. One of our 
motivation for the present work is to define an exact setup to solve this problem for one of the most paradig- 
matic integrable quantum field theory, namely the sine-Gordon model. As a first step we will consider the 
discretized version of the sine-Gordon field theory on a finite lattice. Moreover we will simplify its dynam- 
ics by considering the case for which the representation space of the exponents of the field and conjugated 
momentum is an arbitrary finite dimensional cyclic representation of the quantum al gebra , namely the case 
where the parameter q = e* is a root of unity. This case is interesting in its own right Ill28l4l30 1 and we also 



believe that the treatment of the full continuum theory could then be reached by taking the needed limits; at 
least, we expect to be able to identify the main ingredients and structures necessary to reach this goal and 
learn enough to extend it to other models on the lattice or in the continuum. 

It is worth to describe schematically the microscopic approach that we intend to follow to solve integrable 
quantum field theories by the complete characterization of their spectrum and of their dynamics. 

The first goal is the solution of the spectral problem, for the lattice and the continuum theories: 

i) Solution of the spectral problem for the integrable lattice regularization by the construction of the eigen- 
states and eigenvalues of the transfer matrix using the SOV method. 

ii) Reformulation of the spectrum in terms of nonlinear integral equations (of thermodynamic Bethe ansatz 
type) and definition of finite volume quantum field theories in the continuum limit. 

iii) Derivation of the S matrix and particle description of the spectrum in the infinite volume (IR) limit. 

iv) Derivation of the renormalization group fixed point conformal spectrum in the UV limit. 

The second goal is the solution of the dynamics along the following steps: 

i) Determinant formulae for the scalar product of states in particular involving transfer matrix eigenstates. 

ii) Reconstruction of the local operators in terms of the quantum separate variables. 

iii) Computation of matrix elements of local operators in the eigenstates basis of the transfer matrix. 

iv) Thermodynamic behavior of the above quantities and computation of the physical correlation functions. 



3 An important role in these studies has be en playe d by the fermionic representations of the characters, as derived for different 
classes of rational conformal field theories in Illl7l4l23ll . 

4 A part for s ome of them, like the components of the stress energy tensor, which can be characterized by physical prescription 
flijl and fl2f3l . 
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In this paper, we develop partly this program for the lattice regularization of the sine-Gordon model while 
the completion of it taking into account the required limits to the continuum theory in finite and infinite 
volume case will be addressed in future publications. Let us remark that the possibility to apply the SOV 
method for the discretized ve rsion of the sine-Gordon model in finite dimensional cyclic representations 
was demonstrated recently in Il28ll . hence opening the question of the computations of matrix elements of 
local fields in a completely controlled (finite dimensional) setting. The first result of the present paper is 
to show that the scalar products of states can be computed in this case as finite dimensional determinants 
involving in particular, for eigenstates of the Hamiltonian, the corresponding eigenvalues of the Baxter Q- 
operator; orthogonality of different eigenstates can be proven directly from these expressions. Further, we 
will show that the lattice discretization of the local fields of the sine-Gordon model can be reconstructed 
explicitly in terms of the separate variables. These two ingredients finally lead to the determination of the 
matrix elements of the exponential of the local fields of the model between arbitrary eigenstates of the 
Hamiltonian. 

This article is organized as follows. In Section 2 we define the sine-Gordon model on a finite lattice in the 
cyclic representations and recall the main ingredients of the SOV method in that context. In Section 3 we 
show how to compute the scalar products of states in the SOV representations. The next section is devoted 
to the reconstruction of the local fields in terms of the separate variables. In Section 5 we use these results 
to compute the form factors of the local fields in terms of finite size determinants. In the last section we 
comment on these results and compare them to the existing literature. 



2 The sine-Gordon model 



We use this section to recall the main results derived in II128L 112911 on the spectrum description of the lattice 
sine-Gordon model. 



2.1 Definitions 
2.1.1 Classical model 

The classical sine-Gordon model can be characterized by the following Hamiltonian density: 

Hsc = (d x (j)) 2 + n 2 + 8tt/i cos 2(3cj) (2.1) 

where the field <p(x, t) is defined for (x, t) G [0, R] x R with periodic boundary conditions 4>{x + R,t) = 
c/)(x, t). The dynamics of the model in the Hamiltonian formalism is defined in terms of variables <p(x, t), 
n(x, t) with the following Poisson brackets: 

{U(x,t),(f>(y,t)} = 2TTS(x-y). (2.2) 

The classical integrability of the sine-Gordon model is assured thanks to the representation of the equation 
of motion by a zerocurvature condition: 

[8 t - V(x, t; X),d x - U(x, t; A)] = 0, (2.3) 
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where, by using the Pauli matrices, we have defined: 

U = kicii cos /3<f> + k2<72 sin {}<f> — k3<73±I, (2.4) 
V = — k2CJi cos j3<j) — ki<72 sin j3<j) — ksa^dx^, (2.5) 

ki = i/3(vr^ 1 / 2 (A-A- 1 ), k2 = i/3(vr / /) 1 / 2 (A + A- 1 ), k 3 = ^. (2.6) 
2.1.2 Quantum lattice regularization 

In order to regularize the ultraviolet divergences that arise in the quantization of the model a lattice dis- 
cretization is introduced. The field variables are discretized according to the standard recipe: 

n = 0(nA), n n = An(nA), (2.7) 

where A = R/N is the lattice spacing. Then, the canonical quantization is defined by imposing that <j) n and 
Il n are self-adjoint operators satisfying the commutation relations: 

[<f> n , U m ] = 2m8 nim . (2.8) 

The quantum lattice regularization of the sine-Gordon models that we use here g oes back to [13, 131, 132 1 



and is related to formulations which have more recently been studied in H 133L 1 134TI . Here, the basic operators 
are the exponentials of the fields and of the conjugate momenta: 

v n = e -Wn ; Un = e ^n„/2_ (29) 

Then, the commutation relations (12.81 ) imply that the couples of unitary operators (u n ,v n ) generate N inde- 
pendent Weyl algebras W n : 

u n vm = /" m v m u„ , q = e~^\ (2.10) 



for any n G {1, N}. In terms of these basic operators the quantum lattice sine-Gordon model can be 
characterized by the following Lax matrix^: 

I /yi_„ I Un(g~ 1/2 v„/c„ +<7 1/2 v~ 1 ft- 1 ) (A„v n - (v n \ n y l )/i \ 

Ln[A) - Kn [ (A„/v„ - v n /A n )/i tftf /2 Vntf+<r 1/2 Vn 1 ) 

where A n = A/^ n for any n G {1, N} with ^ n and K n parameters of the model. 
The monodromy matrix M(A) is defined by: 



5 The integrability of the time evolution in the lattice sine-Gordon model is known II 1 3 lL 1 1 35M1 37TI and the most convenient 
way to formulate it uses the Baxter Q-operators; these operators have been constructed for the closely related chiral Potts model in 
113811 . see for a review and rederivation of these points in the lattice sine-Gordon model. 

Here, we make our analysis for the lattice sine-Gordon model leaving the n n and as free parameters. However, the sine- 
Gordon model in the continuum limit is reproduced taking suitable limits on these parameters by relating them to the mass /i and 
the radius R of the compactified space direction of the model. 
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and satisfies the quadratic commutation relations: 

R(X/p) (M(A) 8) 1) (1 <g> MO)) = (1 <g> M(») (M(A) ® l)i2(A//x) 
with the 6-vertex i?-matrix: 



(2.13) 



R(X) 



( qX 



\ 



A - A" 1 g-g- 1 
g-g" 1 A -A" 1 



qX- q X A 1 y 



(2.14) 



The elements of M (A) are the generators of the so-called Yang-Baxter algebra and the commutation relations 
(12.131 ) imply the mutual commutativity of the elements of the following one-parameter family of operators: 



T(A)^tr c2 M(A), 



(2.15) 



the so-called transfer matrix. In the case of a lattice with N even quantum sites, we can also introduce the 
operator: 



e = II» 



(-1)" 



(2.16) 



n=l 



which plays the role of a grading operator in the Yang-Baxter algebra: 



Proposition 6 of 1 128] For the even chain, the charge O commutes with the transfer matrix and satisfies the 
following commutation relations with the generators of the Yang-Baxter algebra: 

9C(A) = c?C(A)0, [A(A),B]=0, (2.17) 
B(A)6 = c?eB(A), [D(A),G]=0. (2.18) 

Moreover, the Q-charge allows to express the asymptotics of the leading generators A(A) and D(A): 



lim A ±N A(A) = ( TT lim A ±N D(A) = ( TT 
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±1 



(2.19) 



and hence the asymptotic behavior of the transfer matrix: 



lim A ±N T(A) 

log A— f=F°° 



n 

\a=l 



±1 



fe + e 



(2.20) 



2.2 Cyclic representations 

In the present article, we will consider representations where both v m and u ra have discrete spectra; in 
particular, we will restrict our attention to the case in which q is a root of unity: 

P 2 = —, P,p' G ^ >0 , (2-21) 
P 
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with p odd and p' even so that q p = 1. The condition (12.211 ) implies that the powers p of the generators u n 
and v n are central elements of each Weyl algebra W n . In this case, we can associate a p-dimensional linear 
space R n to any site n of the chain and we can define on it a cyclic representation of W n as follows: 

v n |&n) = v n q kn \k n ), u n \k n ) = u n \k n - 1), Mk n € {0, ...,p — 1}, (2.22) 
with the cyclicity condition, 

\k n +p) = \k n ). (2.23) 

The vectors \k n ) define a v n -eigenbasis of the local space R n and the parameters v n and u n characterize the 
central elements of the W n -representation: 

< = u£ = <. (2.24) 

Here, we require the unitarity of the generators u n and v n and the orthonormality of the elements of the 
v n -eigenbasis w.r.t. the scalar product introduced in R„. 

Let L ra be the linear space dual of R n and let (k n \ be the elements of the dual basis defined by: 

{kn\k' n ) = {\k n ),\k' n )) = 5 knK Vk n ,k' n e{0,...,p-1}. (2.25) 

From the unitarity of the generators u n and v n , the covectors (k n \ define a v n -eigenbasis in the dual space 
L n and the following left representation of Weyl algebra W n : 

(k n \y n = v n q kn (k n \, (k n \u n = u n (k n + 1\, \/k n € {0, ...,p — 1}, (2.26) 

with cyclicity condition: 

(k n \ = (k n +p\- (2.27) 

In the left and right linear spaces: 

£ N = ^ =1 L n; ^ N = ®" =1 R n , (2.28) 

the representations of the Weyl algebras W n induce cyclic left and right representations of dimension p N of 
the monodromy matrix elements, i.e. cyclic representations of the Yang-Baxter algebra. Such representa- 
tions are characterized by the 4N parameters k = (ki, . . . , k^), £ = . . . ,£n)> u = (""lj • • • > n i\i) an( l 
v = (ui, . . . , un). The unitarity of the Weyl algebra generators u n and v n implies that the parameters v 
and u are pure phases and the following Hermitian conjugation properties of the generators of Yang-Baxter 
algebra hold: 

Lemma 1 of [ 129] If the parameters of the representation and £^ are real for any n = 1, . . . , N and 

satisfy the constrains 

e = —{n n in)l { K nin) * J uniform along the chain, (2.29) 

then it holds: 

{ ' - { Ct(A) Dt(A) ) \ B(eA*) A(A*) J ' ( } 
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which, in particular, implies the self-adjointness of the transfer matrix T(A)/or real A. 
Let us define the average value O of the elements of the monodromy matrix M(A) as 

O(A) = ]JO(q k \), A = X p , (2.31) 

fc=i 

where can be A, B, C or D and we have to remark that the commutativity of each family of operators 
A(A), B(A), C(A) and D(A) implies that the corresponding average values are functions of A, so that B{A), 
C(A) are Laurent polynomials of degree [N] while A(A), T>(A) are Laurent polynomials of degree N in A, 
where we are using the notations: 

r,.i h I 1 for N even, 

[N] = N-e N , N = N + e N -l, e N = { (2.32) 



for N odd. 



Proposition 2.1. 



M(A) 



a) The average values A(A), B(A), C(A), T>(A) of the monodromy matrix elements are central elements. 
Furthermore, they satisfy the following relations: 

{A{A)Y = V{A*), (B(A))* = eC(A*), (2.33) 

under complex conjugation in the case of self-adjoint representations. 

b) Let 

( A(A) B(A) \ 
^ C(A) V(A) ) 

be the 2x2 matrix whose elements are the average values of the elements of the monodromy matrix 
M(A), then we have, 

M(A) = £ N (A)£ N _ 1 (A) ... A (A), 

where 

r ( \\ = ( <? /2 <«* + v n P ) nU^n/a - e n /Av p n )/iP \ 

^ lAJ ~ [ &(A/vZ& ~ ttv P n/A)/iP qP/V(4 n Vn P + V P n ) )' ^ 

is the 2x2 matrix whose elements are the average values of the elements of the Lax matrix L n (A). 



A similar statement was first proven in II139[1 . In the present paper, we will be mainly restricted to the case: 

u n = 1, v n = 1 for n = 1, . . . , N. (2.35) 

In these representations it holds: 

Lemma 2.1. The power 2p of the zeros of B{\) are real numbers with possible complex conjugate couples. 

Proof. The previous proposition and the equality 

C(A) = B(A), (2.36) 
which holds for the sine-Gordon representations with = v„ = 1, as proven in imply that 

(£>(A))* = eB(A*), (2.37) 
and so the statement of the lemma. □ 
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2.3 SOV-representations of the Yang-Baxter algebra 



According to Sklyanin's method [88-90D, a separation of variables (SOV) representation for the spectral 
problem of the transfer matrix T(A) is defined as a representation where the commutative family of operators 

B(A) is diagonal. In il28TI . the following theorem has been shown: 

I — n 

Theorem 1 of [128] For almost all the values of the parameters k and £ of the representation, there exists a 
SOV representation for the sine-Gordon model, i.e. B(A) is diagonalizable and with simple spectrum. 

The proof of this has been given by a recursive construction of the left cyclic SOV-representations for the 
sine-Gordon model. Let us recall here the left SOV-representations of the generators of the Yang-Baxter 
algebra. The Proposition 12. 1 I fixes the average values of B(A): 

/ n \P [N] 

\n=l % J o=l 

in terms of the parameters of the representations. Note that the simplicity of the spectrum of B(A) is 
equivalent to the requirement Z a ^ Zy, for any o / 6 £ {1, . . . , N - e^}. Moreover, the reality condition of 
the polynomial £>(A), proven in Lemma I2TT1 implies that we can chose a N-tupla {r/[°\ ...,77^} of p-roots 
of {Zi, which satisfy the requirements: 

^°)) 2 G^or ^y^K^ 36G{l,...,N}\a: ( ?? (°)) 2 = ((^ 0) )y. (2.39) 
Let {r][ kl \ r]^ N ^\ G £n be the generic element of a basis of left eigenstates of B(A): 

/ N \ [N] 

<k|B(A) = b k (A)(k| with b k (A) = II - ^ n)SN II( A /^ a) " vi ka) /V, (2-40) 

\n=l 1 J o=l 

where we have used the notation 

(k| = (^),...,^ fcN) |, ^)=gSi 0) VaG{l 5 ...,N}, k = (fc 1 ,...,fc N )eZj l . (2.41) 
Then the left action of the other Yang-Baxter generators reads: 

[N] ^ ) 

(k|A(A) = e N ^(-^-(k|T N - ^(k|T+) + £ II tl I, °tf-')(k|T-, (2.42) 



■n (kb) 
( A _ % ■ 

V (k b ) X 



<"|D(A) = e N ^(-^<k|T+ - 2t2 (k |T N) + ^ JT <tf-')(k|T+, (2.43) 



where ry^.A an d f?k.D w& defined by: 



% ° Va 



n N f 



rriv-i ^ 
lln=l W 
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and the shift operators are defined by: 

( V ^\..., V t\...,rj^\Jt = U kl \..,vt ±l \--.,vl kN) \- (2-45) 
The operator family C(A) is uniquel)0 defined by the quantum determinant relation: 

det q M(A) = A(A)D(g~ 1 A) - B(\)C(q~ l \), (2.46) 
where det q M(A) is a central element of the Yang-Baxter algebra (12.131 ) which reads: 

N 

det q M(A) = Yl K 2 n (X/n ni+ - fi ni+ /\)(\/Hn - ~ Mn-A), (2.47) 
n=l 

where [i riy ± = iK^q 1 ^ 2 ^. In the representations which satisfy (12.351 ) the coefficients of the SOV- 
representations can be fixed by introducing the following Laurent polynomials: 

N 

«(A) = n?( 1+ 2<T 1/2 W(l + iq- 1/2 \ n /K n ), d(X) = q N a(-Xq), (2.48) 

n=l 

Note that these coefficients are related to the quantum determinant by: 

det q M(A) = a{X)d(X/q). (2.49) 

Note that for the choice of the parameters {k ti } G iR N and {£ n } G R N the numbers K= Yin=i K n/i an d 
[F n ± are real. 



2.4 SOV-characterization of the T-spectrum 

Let us denote by Sj the set of the eigenvalue functions t(X) of the transfer matrix T(A). From the definitions 
(12.111) and (12.151) . Sj is a subset of M[A 2 , A -2 ]^ where x _1 ]m denotes the linear space over the field 
R of the real Laurent polynomials of degree M in the variable x: 

(f{x))* = f(x*) Vx G C with f(x) G R[x,x _1 ] M . (2.50) 



Note that in the case of N even, the 0-charge naturally induces the grading Sj = U^Lo^ 2 wnere: 



t(\) G S T : lim A 

log A— *-=poo 



±N 



t(A) 



n 

va=l 



±1 



(q k + q~ 



(2.51) 



This simply follows from the asymptotics of T(A) and from its commutativity with ©. 



7 Note that the operator B(A) is invertible except for A which coincides with one of its zeros, so in general C(A) is defined by 



( 12.46b just inverting B(A). This is enough to fix in an unique way C(A), as it is a Laurent polynomial of degree [N] in A. 



The centrality of the quantum determinant in the Yang-Baxter algebra was first discovered in [140], see also [141] for an 
historical note. 
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In the SOV-representations the spectral problem for T(A) is reduced to the following discrete system of 
Baxter-like equations in the wave-function ^(k) = (k| t ) of a T-eigenstate 1 1 ): 



t(4 k ^ t (k) = a(r,^)%(J~(k)) + d(r?( fc "))^(T+(k)) Vr G {1, [N]}. 



(2.52) 



Here, we have denoted by T^(k) = (k±, . . . , k r ± 1, . . . , &n) and a(r]r ) and d(rjr ) the coefficients of 
the SOV-representation as defined in (I2.48l l. In the case of N even, from the asymptotics of T(A) given in 
(12.20b . we have to add to the system (12.521 ) the following equations in the variable rj^ N ^: 



*t,± fc (T+(k)) = <? Tfc *i, ±fc (k) 



(2.53) 



for the wave function \P t ±&(k) = (k| i±fc), where | t±t l ) is a simultaneous eigenstate of T(A) and corre- 
sponding to i(A) G Sj and 0-eigenvalue q ±k with G {0, (p — l)/2}. 



In the paper [129], a complete characterization of the T-spectrum (eigenvalues and eigenstates) has been 
given in terms of a certain class of polynomial solutions of a given functional equation. Let us recall here 
these results; to this aim let us introduce the one-parameter family D(X) of p x p matrices: 



/ *(A) 
-a(q\) 





D(X) 



-d(X) 
%A) -d(qX) 







-o(A) \ 




a(g p - 2 A) t(g p - 2 A) 



V-d^-U) 



p-2 ' 





-d(qP- 2 X) 











(2.54) 



»(g p_1 A) t(^ -1 A) / 



where for now t(X) is just a real even Laurent polynomial of degree N in A. Then the determinant of the 
matrix D(X) is an even Laurent polynomial of maximal degree N in A = X p and we have the following 
complete characterization of the transfer matrix spectrum: 

Proposition 2.2. The set Ej coincides with the set of all the t(X) G M[A 2 , A~ 2 ]^ 2 solutions of the functional 



equation: 



det D(A) =0, VA G C. 
v 



(2.55) 



Moreover, for N odd the spectrum ofT(X) is simple and the eigenstate \ t) corresponding to t(X) G £y is 
characterized by: 

N 

* t (k)^U k i\...,r,^\t) = llQ t (4 kr) ), 



(2.56) 



r=l 



while for N even the simultaneous spectrum of T(A) <2«<f simple and the eigenstate \ t±k) corresponding 
to t(X) G Sj a«<i Q-eigenvalue q ±k with k G {0, (p — l)/2} characterized by: 



N-l 



(2.57) 



r=l 
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Here, 

(p-l)N-(6t+ot) 

Qt(\) = \ at [J ( A ^" A )' 0<a t < P -l, 0<6t<(p-l)N, (2.58) 
ft=i 

a t = 0, 6 t = modp, /or N odd (2.59) 
at = ±& modp, 6 t = ±/c modp, for N even and t(X) G Sj, (2.60) 

f/ze unique (up to quasi-constants) real polynomial solution of the Baxter functional equation: 



t(X)Q t (X) = a(X)Q t (Xg- L ) + d(X)Q t (Xq) VA € C, 



(2.61) 



corresponding to the giventiX), which has been constructed in terms of the cofactors of the matrix D(A), 
in Theorems 2 and 3 of Kl2% . 



3 Decomposition of the identity in the SOV-basis 



The main goal of this section is to obtain the decomposition of the identity in the SOV basis. This is an 
important step towards the decomposition of any correlation function in terms of matrix elements of local 
operators in this basis. In the previous sections we recalled how to define a left SOV basis together with the 
left action of the Yang-Baxter operators on it. To achieve the identity decomposition we need to define a 
corresponding right SOV basis. 

Traditionally, the right states are constructed from the left ones by means of hermitian conjugation. In the 
case of hermitian B operators, this procedure would lead to right states that are still eigenstates of B by 
means of its right action. However, from (I2.30l l. we know that the hermitian conjugate of the B operator 
is the operator C not equal to B. Hence such a procedure (hermitian conjugation of the left SOV basis) 
would lead to a right SOV basis with respect to the C operator, namely it would be an eigenstate basis for 
C. Although such a route could eventually lead to an interesting decomposition of the identity, it would 
result in a very complicated structure of the scalar product of states (in particular, left and right states 
would have no simple orthogonality properties); moreover, we would not be able to obtain, at least by 
obvious means, the coefficients in such an identity decomposition in terms of simple determinants. As a 
consequence, following this path, the structure of the matrix elements of the local operators would hardly 
be given in terms of determinants as these are usually obtained as rather simple and localized deformations 
of the corresponding scalar product formula. 

Therefore, our strategy will be here quite different : we will construct the right SOV basis such that it will 
be an eigenstate basis for the right action of the B operators. Hence it will not be obtained as the hermitian 
conjugate of the left SOV basis; however, due to the simplicity of the B-spectrum, left states will admit 
natural orthogonality properties with respect to right states. Hence, this procedure will lead to a decompo- 
sition of the identity as a single sum over the SOV basis, with coefficients that are not necessarily positive 
numbers but that are computable by means of rather simple and universal determinant formula. To show 
this, we will have to compute the action of left B-eigenstates (covectors) on right B-eigenstates (vectors); 
up to an overall constant these actions are completely fixed by the left and right SOV-representations of the 
Yang-Baxter algebras when the gauges in the SOV-representations are chosen. 
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Let us first define the right SOV-representation with respect to the B operator by the following actions: 



B(A)|k) = |k)b k (A), (3.1) 







A 






A 






_(feo) 
'la 






Jka) 




'la 



A(A)|k) = e N (T+|k)-^- - T S |k)^)5iW + ET + |k) jj O a( „(M), (3.2) 



[ N ] r * ) 

D(A)|k) = eN (T N |k)-^- -Tjlk)^)!^ +^T-|k> FT 4 ^'). 0-3) 

on the generic right B-eigenstate |k) = \rj[ kl \ ...,rj^ N ^) G TZ^. Here, the coefficients a(rj a ) and d(rj a ) of 
the representation are fixed up to the gauge by the condition: 

det q M(A) = d(X)a(X/q); (3.4) 

while C(A) is uniquely defined by the quantum determinant relation (|2.47p . 



3.1 Action of left B-eigenstates on right B-eigenstates 

It is worth remarking that both for the right and left SOV-representations we are explicitly asking that 
the coefficients of the representations of A(A) and D(A) depend only on the zeros of B(A) on which the 
corresponding shift operator acts; i.e. they are separated w.r.t. the zeros of B(A). Naturally such requirement 
is compatible with the Yang-Baxter algebra and the quantum determinant relation but it is not implied by 
them. The main point that we are going to prove here is that this separation requirement completely fixes 
the actions of the generic covector in the left SOV-basis on the generic vector in the right SOV-basis. It may 
be helpful to explain this last statement in terms of the properties of the matrices which define the change of 
basis to the SOV-basis. Let us define the following isomorphism: 

N 

x : (hi, h N ) G {1, ...,p} N -> j =>c(hi,..., h N ) = hi +Y^P {a - 1 \K - 1) G {1, ...,p N }, (3.5) 

a=2 

then we can write: 

(y 3 \ = te|trW = u^ {xi \ and | y .) = u(R)\ Xj) = £ ^)| X . ); ( 3. 6) 

i=l i=l 

where we have used the notations: 

(y^^,...,^! and \ yj ) = \r,[ h l) , r,^), for j = x (h u h N ) , (3.7) 
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to represent, respectively, the states of the left and right SOV-basis and: 

(xj\ = ®n=i{h n \ and \xj) = (g>^ =1 \h n ), for j = x(hi, h N ) , (3.8) 

to represent, respectively, the states of the left and right original v n -orthonormal basis. Here, and 
are the p N x p N matrices for which it holds: 

U {L) B(X) = A B (A)E/ (i) , B(A)[/ (i?) = U {R) A B (X), (3.9) 

where Ab(A) is a diagonal p N x p N matrix. The diagonalizability and simplicity of the B-spectrum imply 
the invertibility of the matrices fTW and and the fact that all the diagonal entries of Ab(A) are Laurent 
polynomials in A with different zeros. Then the following proposition holds: 

Proposition 3.1. Right and left SOV-basis are right and left B-eigenbasis such that the p N x p N matrix: 

M = U {L) U {R) (3.10) 

is diagonal and characterized by: 

M---(ylv -)-(r> (hl) v {hN) \v {hl) JMn _ C " Uali I~fc=i a(^ a) )/atf a " 1) ) 

lll<6<a<[N]Wa 1% ~ % /Va ) 

where we have denoted: 

j = x(hi,...,h n ) (3.12) 
and Cn is a constant characteristic of the chosen representations. 

Proof. The fact that the matrix M is diagonal is a trivial consequence of the orthogonality of left and right 
eigenstates corresponding to different eigenvalue of B(A) and of the simplicity of the B-spectrum. Indeed, 
defined i = x(fci, ku) and j = x(h\, ...,hu), the following identities hold: 

b k (A)My = (^.....^^(A)^,...,^^, (3.13) 
bh(A)M« = (^.....^IBCA)!^,...,^), (3.14) 

where the identity (|3.13p follows acting with B(A) on the left while the identity (|3.14| follows acting with 
B(A) on the right. Then in the case (k\, kn) / (hi, /in) the identities (|3.13p and (|3.14p imply that 
Mij = for any i± j G {l,...,p N }. 

Now, let us compute the matrix element a = (rj^ , ...,rja, 77^ \^(Va)\v[ \ ■•■,Va > \ ■■■> r l^)> wnere 
a G {1, [N]}. Then using the left action of the operator A(fj a ) we get: 

= a^X^, -,»A 0) |t|J 0) , ..-,11^, -..,^ 0) > (3-15) 

while using the right action of the operator A(?7 a ) and the orthogonality of right and left B-eigenstates 
corresponding to different eigenvalues we get: 

I N 1 1 (i) / (o) (o) , (lh 

b^a,b=liVa /% —f] h /T)a ) 



17 



and so: 



(r^r^ M tf>/^-^/^) _ (31?) 

Then by applying (13.171) /i a times, we get: 

(r ? (Q) 1 ...,^),...,^|r ? W,...,^,...,^) _^ q(^)) ™ (^V^-^/^) 

(3.18) 

Now, let us consider explicitly the case of even N. In this case we still have to define the recurrence for 
a = N. We compute the matrix element: 

N (A) ee (^.....^.^lACA)!^,...,^,^), (3.19) 
acting with A(A) on the right we get: 

while acting on the left we get: 

77m I7v a 



»7n ^k.A 

which simply implies 



(3.21) 



(3.22) 



The previous formula implies, for both N even and odd: 



/„(0) JO) JO), (0) (0) (0)v 11 11 3r _(fea-l)^ 11 ._(*.) " 

(3.23) 

To prove it let us observe that (13.181 ) coincides with (13.231 ) for = ... = = 0. Now in the case 
/13 = ... = /iM = with hi ^ and /12 7^ 0, we get (13.231) by the following factorization: 

(v^\4 h2 \--A hN) H hl) A h2 \--A hN) ) 

Aft aj«th n (W-W) 

(3.24) 
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and so on for the generic case. Finally, from (13.231 the statement of the proposition follows being by 
definition: 

Ma _ (^\vi h2 \..A hN) H hl \^\..A hN) ) 

□ 

Remark 2. The diagonal matrix M is mainly fixed by the requirement that the left and right representations 
have a separate form in the zeros of B(A). Indeed, this fixes completely the denominator in all the entries of 
M. Moreover, the constant Cm is a function of the representation only via the central elements {Z\, Z^). 
In the following, we will fix: 

CN = rfV /2 ) eN , (3-26) 

this choice just amounts to an overall renormalization of the states. Finally, the products of a(\)/a(\q~ 1 ) in 
the numerator of each matrix element is fixed from the choice of the gauge done in the SOV-representations. 
Note that we are always free to take the following gauge: 

aiXq- 1 ) = o(A), (3.27) 

for which the numerator in (13.1 11 1 is Cn . 



3.2 SOV-decomposition of the identity 

The previous results allow to write the following spectral decomposition of the identity I: 



i = I>|yi><yil, (3-28) 



i=i 

in terms of the left and right SOV-basis. Here, 



Hi = -. — j — r-, (3.29) 

is the required analogue of Sklyanin's measur^, which is discrete for the cyclic representations of the sine- 
Gordon model. Explicitly, the SOV-decomposition of the identity reads: 



n rr[ N l ( ( h 

hi,...,A N =l l<fe<a<[Nl °N llb=l UbKVb 



i= e n (^ a) f-^ b) ) 2 ) - ( ;;;; ,7?N (3-30) 



where the separate functions cJbivi ) reproduce the numerator in (13.11) : 

fNl-1 

k b =l 



- 6 tf 6) )-(^ ) ) [N1 " 1 n^ ) )/a(^- 1) ) (3.31) 



9 Note that here it cannot be called a measure as the B(A) operator is not self-adjoint, meaning that in general the fa are not 
positive real numbers. However the above procedure indeed defines a proper decomposition of the identity operator with simple 
and computable coefficients. 
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and they are gauge dependent parameters. 

Remark 3. Sklyanin's measure0has been first introduced by Sklyanin in his article on quantum Toda chain 
|88[]. There, it has been derived as a consequence of the self-adjointness of the transfer matrix w.r.t. the scalar 
product. In particular, the Hermitian properties of the operator zeros and their conjugate shift operators have 
been fixed to assure the self-adjointness of the transfer matrix. In the similar but more involved sinh-Gordon 
model niOOll . the problem related to the uniqueness of the definition of this measure has been analyzed. 
There, it has been proven that the measure is in fact uniquely determined once the positive self-adjointness 
of the generators A(A) and D(A) is required. In the compact case of the sine-Gordon model the method 
used here insures that the analogue of Sklyanin's measure is uniquely determined up to an overall constant 
and the choice of gauge, as d iscus sed in the previous remark. Let us mention that the S klyanin's measure 
has already been derived in 1 143 1 for cyclic representations of the related T2-mode0 Il47 - 149 1. There 
the measure has been obtained by a different approach, i.e. by a recursive construction which needs to go 
through the recursion in the construction of left and right S OV-basis. It is interesting to remark that in our 
purely algebraic derivation we skip these model dependent features so that the approach we used is suitable 
for general compact SOV-representations of the 6- vertex Yang-Baxter algebra. 



3.3 SOV-representation of left and right T-eigenstates 

Up to an overall normalization, the SOV-decomposition of the identity and the SOV-characterization of the 
transfer matrix spectrum lead to the following representations of the right eigenstate of the transfer matrix 
T(A): 

\t±k) =^2^i{ji\t±k)\ji) 
i 

P / n Tkh N \^ [N] | (hi) (M\ 

e (v) n<2rf j > n (tf-') 2 -^') 2 ) 1 ^]'^; X, ■ < 3 - 32 > 



i=l 



corresponding to the eigenvalue t(X) € Ey, where the index k € {0, (p — l)/2} appears only for N 
even and indicates that \t±f.) is also a 0-eigenstate with eigenvalue q ±k . Here Qt{X) is the solution of 
the Baxter equation (13.341 ) defined in Proposition 12.21 In a similar way we can prove that, up to an overall 
normalization, the left T-eigenstate has the following SOV-representation: 



(<±*l = $^M»<*±fc|y»)(ytl 

P /_±M» N \ "N [N] / (hi) (An) I 

= E (V) n<3.("i"-') n tf-'-^'l ^rrrW - (3 - 33 » 

h u ...,h N =l V y 7 a=l l<6<a<[N] lh=l u b\nb ) 

where Qt(X) is the unique (up to quasi-constants) polynomial solution of the Baxter functional equation: 

t(\)Q t (X) = (J(A)Qt(Ag- 1 ) + a(X)Q t (\q). (3.34) 



— 

See also 1 97] for further discussions on the measure. 
"The SOV analysis for these representations has been first developed in Il44ll . 
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Remark 4. In the gauge fixed by (13.27I ). this Baxter equation reads: 

i(A)Q t (-A) = q N a(X)Q t (-Xq- 1 ) + q~ N d(\)Q t (-\q), (3.35) 
and so, up to quasi-constants, we have: 

Qt(X) = X XN Qt(~X), X n = N modp, 0< X n <P-1- (3.36) 

4 Decomposition of the identity in the T-eigenbasis 

Here we use the results of the previous section to compute the action of covectors on vectors which in the 
left and right SOV-basis have a separate form similar to that of the transfer matrix eigenstates. 

4.1 Action of left separate states on right separate states 

Let us start giving the following definition of a left (a^ | and a right separate states; they are respectively 
a covector and a vector which have the following separate form in terms of the SOV-decomposition of the 
identity: 

P /Jfc/l N \ e N / \2 

w= £ (hp) n-hf-') n ((■£■>) - 

/il,...,ft N =l 7 o=l l<6<a<[N] 

i«= e (V) n n (("?"') 

/ii,...,h N =l V y 7 0=1 l<b<a<[N] 

where the index k appears only for N even. It is easy to see that such states generate the whole space of 
states of the model (in particular the T-eigenbasis is just of this type). The interest toward these kind of 
states is due to the following: 

Proposition 4.1. Let (a^| and \/3k') two separate states of the form (|4.ip and ([4. 21) . respectively, then it 
holds: 

MM = ^,def *H ■** - W*"*± ^P ^* (4-3) 

lN| h=i ^o(*?o ') 

Proof. From the SOV-decomposition, we have: 

\h N =l F J hi,...,h [N] =l 0=1 w a(,f?a ) 

where V{x\, xn) = ni<6<a<[N]( x o ~~ x &) * s tne Vandermonde determinant. From this formula by using 
the multilinearity of the determinant w.r.t. the rows we prove the proposition. □ 



(hi) 



_(h N )| 
i 'In 



rrt N ) /• 

lUJi^Wb 



(h b 
% 



b ) 
(hi) J^ti)\ 



)\\\rir,..,vr) 



(hi). 



(4.1) 



(4.2) 
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In Section I2T21 we have associated to the linear space 7^n tne structure of an Hilbert space by introducing a 
scalar product. Then it is clear that the above determinant formula represents also the formula for the scalar 
product of the two vectors ((afc|)^ and \/3h) in7£|\|. It is worth pointing out that the vector ((afel) G T^n is a 
separate state in the right C-eigenbasis, as it simply follows from the hermitian conjugation properties of the 
Yang-Baxter generators reported in Section [2721 Then these results can be considered as the SOV analogue 



of the scalar product formulae [27 U142H computed in the framework of the algebraic Bethe ansatz. However, 
we want to stress that the determinant formulae obtained here are not restricted to the case in which one of 
the two states is an eigenstate of the transfer matrix, on the contrary to what happens for the scalar product 
formulae in the framework of the algebraic Bethe ansatz. Finally, we can prove directly from these formula 
the following: 

Corollary 4.1. Transfer matrix eigenstates corresponding to different eigenvalues are orthogonal states. 

Proof. Let us denote with \tk) and \t' h ) two eigenstates of T(A) with eigenvalues i&(A) and t' h (X) for N odd 
and with eigenvalues q h and q h for N even. To prove the corollary, we have to prove that: 



detail = with <f S ,<»> ^ y^ , (4.5, 

with h = k for N even. To prove (I3.27I ). it is enough to show the existence of a non-zero vector V(*>* ) such 
that: 



E*SM* =° VaG{l,...,[N]}. (4.6) 



v (/./'. v !/./', 
6=1 

For simplicity, we construct this vector in the gauge (13.271 where it results: 



Let us recall that the eigenvalues of the transfer matrix are even Laurent polynomials of degree N of the 
form: 



t*(A) = eNln^ltf + fV + A-VE^* 1 ^ ( 4 - 8 ) 
\a=l Z / 6=1 

/ N ±1 \ [N] 

4( A ) = e N (/ + ^)(A N + A- N ) + ^c b A-[ N ]- 1 + 2& , (4.9) 

\a=l * / 6=1 

so if we define: 

Vf*' 5 = 4 - Q, V6 G {1, [N]}, (4.10) 

it results: 

[N] p 

E^i'fvf ' } =E^'(^ C) )^(^ C) )(^(^ C) ) (4- 11 ) 

6=1 c=l 



±1 
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We can use now the Baxter equations (12.611 ) and (I3.34I ). with the chosen gauge (13.27I ). to rewrite: 

QAvi ha) )Qt(vi ha) )(t'(vi ha) ) -t(v ( a a) )) = a(vi ha) )QAvt- 1] )Qt(rii ha) ) 

+ rf(^ ) )Qt'(r/i /la+1) )Qi(^ ) ) 

- a(^ a+1) )^(^ a) )Q t (^ a+1) ), (4-12) 
and by substituting it in (14.111 ) we get (14.6b . □ 

4.2 Decomposition of the identity in the T-eigenbasis 

Let us remark that the diagonalizability and simplicity of the transfer matrix spectrum implies the following 
decomposition of the identity in the left and right T-eigenbasis: 

'-"if E Hf. 

where 

(tk\h) = det ||*2J^|| with = (V^) 2 ^ £ Qt ^)Q^ c) ) g2( 6-i) C) (4 _ 14) 

c=l ^ai^a ) 

is the action of the covector on the vector \tk) as defined in Section [331 It is worth to note that in the 
representations which defines a normal transfer matrix T(A), the generic covector (tfe|= (\tk))K dual to the 
right T-eigenstate \tk), is itself a left eigenstate of T(A) which taking into account the simplicity of the 
T-spectrum has to satisfy the following identity (ifc|= at k (tk\, where (tk \ is the left T-eigenstate defined in 
(13.331) . Of course, in these representations the following identities hold: 

\tk){t k \ = (4 15) 

(tk\t k ) \\\t k )\\ 2 

where \\\tk)\\ is the positive norm of the eigenvector \tk) in the Hilbert space TZ^. The above discussion 
implies the relevance of computing explicitly the norm of the transfer matrix eigenstates (13.321 ) as it allows 
to fix the relative normalization at k thanks to the identity at k = |||ifc)|| 2 /{tk\tk) and then it allows to take 
these left and right states as the one being the exact dual of the other. This interesting issue is currently 
under analysis. 



5 SOV-representation of local operators 

The determination of the scalar product formulae, presented in the previous section, is the first main step to 
compute matrix elements of local operators. The second one is to get the reconstruction of local operators in 
terms of the generators of the Yang-Baxter algebra, i.e. to invert the map which from the local operators in 
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the Lax matrices leads to the monodromy matrix elements. Indeed, the solution of such an inverse problem 
allows to compute the action of local operators on the eigenstates of the transfer matrix. Together with the 
scalar product formulae it leads to the determination of the matrix elements of local operators. 

The first reconstruction of local operators has been achieved in 02711 . in the case of the XXZ spin 1/2 chain. 



In 129], it has been extended to fundamental lattice models, i.e. those with isomorphic auxiliary and local 
quantum space, for which the monodromy matrix becomes the permutation operator at a special value of 
the spectral parameter. The reconstruction also applies to non-fundamental lattice models, as it was shown 



in 12911 for the higher spin XXX chains by using the fusion procedure [ 150]. In the case of the sine-Gordon 



model this type of reconstruction is still missing and the only known results are those given by T. Oota 



based on the use of quantum projectors 11511. However, it is worth recalling that Oota's results only lead to 
the reconstruction of some local operators of the lattice sine-Gordon model. In this section, we will show 
how to obtain all the local operators of the sine-Gordon model for the cyclic representations which occur at 
rational values of the coupling constant /3 2 . 

5.1 Oota's reconstruction of a class of local operators 



Here we recall some of the results of Oota IU51I1 which lead to the reconstruction of a certain class of local 
operators in the sine-Gordon model. 

The Lax operator L n (A) has the following factorization in terms of quantum projectors: 
i („ \ P n - ~ | u^ /2 (v nKn + v" 1 ^ 1 ) \/ 1/2 _ 1/2 \ 

y U„ (V n /i n x + V n K n ) J \ / 
( 1/2 \ 

K(Pn-) = Pn-Qn- = K n f j ( (v n K n + V^K" 1 ) uj/ 2 (v n K^ 1 + V' 1 K n ) U,^ 2 J , (5.2) 

when computed in the zeros fj, Ui ± of the quantum determinant; such factorization properties are at the basis 
of the following Oota's reconstruction. 

Proposition 5.1. The local operators u n and ao,n = {{Q ly n + K n)/ {l~ ly n K n + l)) u n 1 admit the re- 
constructions: 

u n = U n B- 1 ( / L in , + )A(/i ni+ )U- 1 = U n D- 1 ( Mn , + )C(^ n , + )U- 1 , (5.3) 

a ,n = U n A- 1 (^ ni _)B(^ ni _)U- 1 = U n C- 1 (^ ni _)D(^_)U^ 1 , (5.4) 
where the shift operator U n brings the quantum sites from 1 to n — 1 to the right end of the chain: 

U n Mi,... iN (A)U- 1 = M„,..., N>ll ..., n _i(A) = U_x(A) • • • k(A)L N (A) • • • L n (A). (5.5) 

It is then clear that the formulae (I5.3l l- (I5.4I ) allow to reconstruct all the powers u^ = 
U n (B~ 1 (^i nj+ )A(^ nj+ )) fc U" 1 of the local operators u n but they do not give a direct reconstruction of 
the local operators v n ; indeed, only rational functions like (q~ 1 v'^ l + k^)/ fg _1 v^/c^ + l) are obtained. 
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Let us make some comments on the shift operators U n . The definition (15.51 ) characterizes the shift operator 
U n up to a constant and implies, by the cyclicity invariance of the trace, their commutativity with the transfer 
matrix T(A). Moreover, in the case of even chain, the shift operators U„ clearly commute also with the 8- 
charge. Then, in the cyclic representations of the sine-Gordon model under consideration, the simplicity of 
the transfer matrix spectrum implies: 

Un\tk)=<P ( n k) fa), (5-6) 

where \tk) is the generic eigenstate of T(A) for odd chain and of (T(A),B) for the even chain. In particular, 
this implies that the shift operators only produce a prefactor in the form factors of loc al op erators which is 
one if left and right eigenstates are dual of each otheJ^i It is worth remarking that ir0 Ii32h . for the special 
case of highest weight representations of the even lattice sine-Gordon model, it has been shown that: 

a=l o=l W/M 

where fi a are zeros of the quantum determinant in these representations. This result is interesting as it shows 
that the shift operators U„ for non-fundamental lattice models, like the sine-Gordon model, are characterized 
by the same type of eigenvalues they have in fundamental lattice models, like the XXZ spin 1/2 chain. 
However, the proof of (I5.71 i presented in 111 3211 is representation dependent, as it is based on the algebraic 
Bethe ansatz representation of the transfer matrix eigenstates. Then, an independent proof is requ ired for 
cyclic representations of the sine-Gordon model and it will be given directly in a future publication 015211 in 
the more general cyclic representations of the 6-vertex Yang-Baxter algebra associated to the T2-model. 



5.2 Inverse problem solution for all local operators 

Here, we show how to complete the reconstruction of local operators by solving the inverse problem for the 
local operators v„ and their powers. The main ingredient used will be the cyclicity of the representations of 
the sine-Gordon model here analyzed. 

Let us define the local operators: 

R — ,, fc ™ ,,1-fc _ ly n + K n K m 

Pk,n = U n a ,nU n - 2fc _! 2 2 : , (5.8) 
1 v n K n "r 1 

then the following proposition holds: 

Proposition 5.2. In the cyclic representations of the sine-Gordon model, the local operators v^ fc admit the 
following reconstruction: 

r -\\k / 2p 2p . h\ p-i 
vf = ( \ Vn 2 Kn ^ ] Y.1~ k{2a ~ l) ^ forke{l,...,p-l}, (5.9) 

P K n \ K n K n ) a=0 

where: 

, l-k' 



;B- 1 ( / u n , + )A( Mn , + )) fc A- 1 ( / i„ i _)B(^ i _)(B- 1 ( / x ri , + )A(^ n>+ )) 



(5.10) 



12 Translational invariance for the limit of the homogeneous chain. 
13 This result is there attributed to V. Korepin, private communications. 
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Proof. In our cyclic representations the local operators u n and v n satisfy the property that u„ and are 
central, i.e. \F n and are just numbers vP n and Vn which characterize our representations. Then the following 
identity holds: 

p-i 

n(? 2i_i «+ i)=i+« ; (5.H) 

j=0 



and so: 

2p 2p , 1 P-l 

£^ ^ : = E(- 1 ) a « a(2 ^ 1)v « a ^- ( 5 " 12 ) 

H n n > a=Q 

The previous formula allows to rewrite the rational function ^ (7l as a finite sum in powers of v^: 

„2p 2(p~l) , 2 , ( 2 _ -2\ ( _i\a a(2k-l) 2a 2a 

Pfc,n — 2p 2p -. ' 

Vn ~T t 



then by taking the discrete Fourier transformation, we get the reconstructions (I5.9l l. plus the sum rules: 

P-l 2p 2(p-l) 2 

Finally, the representation (15.101 ) for the /? aj „ are trivially derived by (I5.3t - d5.4l i. □ 
Note that thanks to the identities = v n h /vn for k = 2h — p odd integer smaller than p, the formulae (15.91 ) 



indeed lead to the reconstruction of all the powers for k 6 {1, ...,p — 1}. Then the previous proposition 
together with Oota's reconstructions leads to the announced complete reconstruction of local operators for 
cyclic representations of the sine-Gordon model. 



5.3 SOV-representations of all local operators 

In order to compute the action of the local operators and on eigenstates of the transfer matrix and 
eventually obtain their form factors, it is necessary to first derive their SOV-representations0. Here, we 
show how these SOV-representations can be obtained from the previously given solutions of the inverse 
problem. In order to simplify the presentation, we introduce here explicitly the operators^ 771 , . . . , 77N > "Hh an d 
?7d defined by the following actions on the generic element (k| = (rj[ kl \ rj^ N ^ | of the left B-eigenbasis: 

(k|r? a = r?i fca) (k|, Va € {1, N}, (k^A = %,A(k| and (k|7? D = r ?k , D (k|, (5.15) 

where the corresponding eigenvalues rj a , ry^A an d ?/k,D w& the complex numbers defined in Section 1231 
The following lemma is important as it solves the combinatorial problem related to the computations of the 
SOV-representations of monomials in the Yang-Baxter generators: 

14 To simply the notations we chose to present the results in this subsection only for the case N odd and for the representations 
withv* = u£ = 1. 

15 To simplify the exposition, we have decided to keep as simple as possible the notations for these operators, we hope that 
nevertheless the difference with the corresponding eigenvalues is clear. 
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Lemma 5.1. The operator (B 1 (A)A(A)) has the following left SOV-representation: 



(B- 1 (A)A(A)) fe = K- fc J2 

a={ai,...,a N }6Ng: 



A: 
a 



n n 

j=l \h=0 



a(VjQ h ) 



{Xq h / Vj - ^/Xqh) 



n n 



i^j,i=l h=ai—aj+l 

acting on the state (rji, ...,r}n\, where: 



(VjQ h /Vi ~ Vi/VjQ h ) 



(5.16) 



k 
a 



U!J=i [ay]! 

Proo/ Here we use the commutation relations 



W , [ k ]\ = [k][k-l]-[l], [a] = q " "~ 



q-q- 



(5.17) 



Tim = q ±5a Vi, 



(5.18) 



and the following left SOV-representation: 



B- 1 (A)A(A) = K- 1 ^ 



_rr t 



(5.19) 



then the lemma holds for k = 1 and we prove it by induction for k > 1. Let us take N integers c^: 

N 



= (5.20) 

a=l 

from which we define the set of integers I = {i G {1,...,N} : a; / 0} and c£ fc) as the coeffi- 
cient of n T «~ ai in the expansion of the k-th power of B _1 (A)A(A). By writing (B _1 (A)A(A)) fe = 
(B- 1 (A)A(A)) fe - 1 B" 1 (A)A(A) and by using the induction hypothesis for the power k - 1 of B _1 (A)A(A), 
we have: 



ael 
N aj—8a,j — l 

n n 

j=l h=0 



k - 1 

a - 6 a 



aim h ) 



(Xq h /r ]j -r ]j /Xq h ) 



N 

n 



»o + l^ 



n 



h r]j/r]i — rji/q ai s °- i h rjj 

1 



iG-f \{a} 

with 5 a = . . . , <5|\i,a)- The first term in the r.h.s. is the coefficient of f\ T 

and the second is the coefficient of T" 1 in B _1 (A)A(A) once the commutations between f\ T i ai+5a 



(5.21) 

n,+ ' W in(B- 1 (A)A(A)) fe - 1 
and 
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the rji have been performed. This can be rewritten as: 
k [k-l]\ ( N N 

C ^ = ttr! 1 n n ( n 

111 11 \j=l h=0 &j,i=l H 13 



a{q \) 



n 



lm - i]i/Q ai hr ij (M h /vj - r n/M h ) 
a mha - i]a/q aa i]i . 



ael iel\{a} 

which leads to our result when we use the relation: 

f a mha - rja/Q aa Vi 



i rnha - i]a/q aa ai m 



(5.22) 



0=1 



.a=l 



(5.23) 



The above formula holds for any n, for any set of numbers rji and for any non-negative integers on, all 
of them being in generic position. This is shown by studying the contour integral and the residues of the 
function: 



9{z) 



n 

i=i 



-2a ir j2 ' 



(5.24) 
□ 



Remark 5. Let us point out that the power p of B 1 ( A) A(A) is a central element of the Yang-Baxter algebra 
and it reads: 

(B- 1 (A)A(A)f = B(A)-U(A), (5.25) 



as it simply follows from the commutations relations: 

B- 1 {qX)A{qX) = A(A)B _1 (A). 



(5.26) 



Then, it is important to verify that the same result follows from the previous lemma for k = p. In order to 
prove it, it is enough to use the following properties of the quantum binomials: 



p 




a 





1 if 3i G {1, N} : at = p6 a ,i Va € {1, N}, 
otherwise, 



(5.27) 



from which the formula (I5.161 l for k = p reads: 



^ (A/Z a - Z a /A) 11 (Z b /Z a - Z a /Z b ) ' 



(5.28) 



and to observe that the r.h.s. of (I5.28l l indeed coincides with B(A) 1 A(A). 

Note that the above lemma gives directly the SOV-representations of the local operators and k n with 
k £ {1, ...,p — 1} when we fix the parameter A to [i n ^ e with e = + and e = — , respectively. Moreover, it 
allows to derive also the SOV-representations of the local operators as it follows: 
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Corollary 5.1. The local operators v^ k with k G {1, ...,p — 1} have the following left SOV-representation: 



1=1 a={a u . ..,a N }€N^. 



p — 1 



N 



,(2fc) 



N 



X II II ( vq h/ vrj / vq h\J y n,(a,a)LL 

i^j,i=lh=a i -(a j +S j , a )+l ' 11 111 m 1 I .7=1 



TT r - 



(5.29) 



v ( 2fc ) = 



(5.30) 



?{2fc) _ (-i)fc(^ + i) ! 



n,(a,a) 



_y - 



-k(2r-l) {K 2 n _ K -2 )a{rjaq -a a) 



+ J .\ 



N r-1 

nn 



(^ + q^ +h h 3 -r ]j /^ + q^ +h ) 
{^n,+q h hj - r]j/p, n ,+q h ) 



(5.31) 



3=1 h=0 

Proof. This is a consequence of the previous lemma and of the identities: 



2 —2 — k 

U- 1 /3 fc , n U n = _,^ w ~_" T L„_ 2 7fc,n + i \ ~ q _, _ 2 , for fc G {1, - 1}, (5.32) 



q k K 2 n - q k K r 



k 2 -2 



where: 



P-i 



fc-i 



7fe , n = B" 1 ^) [J B( /Un , + ^)A- 1 ( /Un ,_)B(/i„,_)B- 1 (/i ri , +(? fe )A(/,„, + ^) B(/in, + ^'), (5.33) 

j=k j=0 



Now by using the relations: 



A- 1 ( A i ni _)B(/j ni _) = (B- 1 ( M „ ) _)A(/j ni _)) 



p-i 



(5.34) 



we get the following representations: 

N 



p-1 

a 



N 

n n 



= E E 

a=la={ai,...,Q N }6NN : L J j=l \ h=0 

ELi a h=P-l 
N «, 

x TT TT (Hn,+q a]+h /Vj ~ Vj I 'fin,+q ai+h ) A -(a J+ ^. a ) 



(t*n,-Q h /Vj -Vj/Vn,-q h ) 



a(rj a q Qa ) 



from which our result follows. 



(5.35) 
□ 
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6 Form factors of local operators 



In the following we will compute matrix elements (form factors) of the forrrj^l: 

(t\O n \t'} (6.1) 

which by definition are the action of a covector (t\ € £n> a left T-eigenstate defined in (13.33b . on the vector 
obtained by the action of the local operator O n on the right T-eigenstate \t') £ IZ^ defined in (I3.32I ). Of 
course, these form factors depend on the normalization of the states (t\ and \t') and then it is worth pointing 
out that nevertheless we can use them to expand m-point functions like: 

— m — ■ <U) 

Indeed, by definition these m-point functions are normalization independent and by using m-1 times the 
decomposition of the identity (|4.13|) . we get the expansions: 

<t\0«,---0«Jt) v (t|0„,|(<")(t'-"|0»„|t)n".'2(' ( '""|0»J' ( ' ,) ) , 63) 

(*> km** ei*>nS<*>l*e»> ' l "' 

where in the r.h.s there are exactly the form factors (|6.ip that we are going to compute in this paper. 



6.1 Form factors of u r 



In this section we use the SOV-representations of local operators to give some examples of completely 
resummed form factors. 

Proposition 6.1. Let (tf-\ and \t' k ,) be two eigenstates of the transfer matrix T(A), then it holds: 



„(**) 



(t k K\t' k ,) = fc' + i) 6N <&(ll<f W)ll)> 



r 



(*fc/) x ' [N] 



(6.4) 



where (fn k> and (fn k " are the eigenvalues of the shift operator U n and \\U^£ is the [N] x [N] matrix: 

(6.5) 



(A) = <£ 1/2 6G {!,..., [N]-!}, 



<rS(A) 



jo) 

'la 



-1 



(0) 



E 



fM-v h QM h) : 



A 



(/l+l)> 



(A/^ +1) -# +1) /A) 

N 



A 



— rrsii 



,-fc' 



(6.6) 



w/iere K= ]ln=i K nA 



16 To simplify the notations in this introduction to Section [6]we are omitting the index k £ {0, (p — l)/2} in the transfer 
matrix eigenstates which are required in the case of even chain. 
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Proof. The right SOV-representation of the operator B (A)A(A) reads 



B- 1 (A)A(A)|k) = ^ (A T + + ^T N ) |k) + £ 



n 



0=1 



Kr ?N N ( A /^g - r] a q/X) Al (ri a /ri b - rj b /r] a ) ' 

(6.7) 

Let us denote with [B _1 (A)A(A)] the second term on the r.h.s. of (I6.7l i. Then, we observe that from the 
SOV-decomposition of the T-eigenstates, we have: 



(t k \[B-\X)A(X)]\t[ 



' V p Jk-l-k')h N \ e N [N] p 

Wl? -EE v (\'ii 

/ 0=1 hi,...,h\fin=l 



(0) 



Oh) 



[N] 

>< n 



vl hb) QAvi hb) )Qt(vi hb) ) 



(6.8) 



and so: 



/ \ e N [N] p 

<t fc |[B- 1 (A)A(A)]|f fc/ )= P^±M £ E fc( 

\ / a=l /n,...,/i N =l ^_ 



- n 

6^a,6=l 



/i a is missing. 

^ b) Qt>{vt b) )Qt(vt b) ) 
Mvl hb) ) 



(We have removed the row a.) 



(_l)[N]+o ^ 



P Qt(vt +1) )QAv { a ha) )Ua ha) 



-l 



a(?7a 



^ Catf a) )K(A/^ +1) -^ +1) /A) 



(6.9) 



bringing the sum over (hi, h a , /i[Mj) inside the Vandermonde determinant V a , we have that the above 
expression is just the expansion of the determinant: 



(t fc |[B- 1 (A)A(A)]|4,) = (* fcifc , +1 )" N det(|| U^piX) 



(6.10) 



where 



Ujtf\\j\ coincides with &%Q 1/2 for b G {1, [N] - 1}, while: 



n (o) 

'/a 



-1 



/ MwW (M-l) /l\ ^ a 



(6.11) 
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Now let us compute the matrix elements: 



/V p n (k-l-k')h N \ e N P 9 



PVn llj=l$j / hi,...,h [N] =l 



2 



^ > ( 6 - 12 ) 

6=1 Mvl ) 



and so: 

By using the fact that [N] — 1 columns are common in the matrix of formula (I6.101 l and in those of (I6.13I ). 
we get our result. Let us remark that the above result holds for any value of A. □ 

Remark 6. It is worth pointing out that the form factors of u n are written in terms of a determinant of a 
matrix whose elements coincide with those of the scalar product, except for the last line which is modified 
by the presence of the local operator. It is then interesting to recall that a similar statement holds for the 
form factors of the local operators in the XXZ spin 1/2 chain. 

6.2 Suitable operator basis for form factor computations 

In this section we introduce an operator basis which can be conveniently used to describe all local opera- 
tors. The interest toward this basis is due to the fact that the form factors of its elements are simple being 
represented by a determinant formula. 

6.2.1 Basis of elementary operators 

Let us introduce the following operators: 

n - B(^ +fc ' 1} )B(^ +fc - 2) ) • • • Btf +1 >)AfoW) , r fn n 
pn^Kip-v U [ S a , b =i(zjz b - z b /z a ) 

where the rj a k ^ are fixed in Section [2731 

Lemma 6.1. The operators Q a ,k satisfy the following properties: 

Oa,kOa,h is non-zero if and only ifh = k — 1, (6.15) 

and 

O a kO a k-~i ■ ■ ■ O a k+i-pO a k- p = — frrj O a k- (6.16) 

Y\ [ Sa, b=1 (Za/Z b -Z b /Z a ) 

Moreover the following commutation relations hold: 

m O a>k = qO atk r, A , [ m ,O ajk ] = [G,O a<k ] = 0, (6.17) 
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and 

n n / % % I 'to )_ n n /-aio\ 

fora^be{l,...,[N}}. 

Proof. The first property follows from B{Z a ) = 0, where £>(A) is the average value of B(A). By the 
definition of O a> k it is clear that: 

/„(*0 Ah) J k ^\n u - a(Va ) )Sh,k , (fei) (fc_i) (An), (619) 

llb^a,b=lV'/a I 'lb % I 'to ) 

so that the second property simply follows. To prove the last property we have to use the Yang-Baxter 
commutation relation: 

(A/// - /*/A)A(A)B(/z) = (A/g/x - W /A)B( / u)A(A) + (g - g" 1 )B(A)A(/i) (6.20) 

we have before to move the A (77^) to the right through all the B(rja ), remarking that only the first term of 
the r.h.s of (16.201 ) survives, and after to move the A(ry„ ) to the left. □ 

Let us introduce now the following monomials which we will call elementary operators: 

J u ai,ki'" u a r ,k r > { ~°- ZL) 
where YIh=i a h < P, k,ki G {0, ...,p - 1}, a, < aj G {1, [N]} fori < j € {1,..., [N]} and: 

°2 = Ca.feCa.fe-l • • • Oa,fc+l-a, with a G {1, ...,p}. (6.22) 
Then the following lemma holds: 

Lemma 6.2. For any n G {1, N}, the set of the elementary operators dressed by the shift operator \J n : 

II f(<Xl,-,<*r) lj-1 ((,21) 

is a fow/s in the space of the local operators at the quantum site n. 

Proof. The space of the local operators in site n is generated by ujjandvjjforfc G {l,...,p— 1}. Proceeding 
as done in Proposition 15.21 we have in particular the possibility to show that an alternative local basis is 
defined by the operators: 

u k n = U n (B- 1 (^n, + )A(/x„, + )) fc U"\ (6.24) 
h,n = U n (B- 1 ( / x n , + )A( / x n , + )) fc B- 1 ( / u n ,_)A( / u n ,_) (B- 1 ( / x„ i+ )A( /Uni+ )) p " 1 " fc IL; 1 (6.25) 

for k G {1, — 1}. Then to prove the lemma we just have to show that the above operators are linear 
combinations of those defined in (16-2 1 b - Note that for A p / Z a with a€{l,...,[N]}, the operator B _1 (A) 
is invertible and by the centrality of the average values we can write: 

B->(A)A(A) _ B(A/-)B(A^)...B(A,)A W ^ 
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Now the operator B(Xq p 1 )B(Xq p 2 ) • • • B(Xq)A(X) is an even Laurent polynomial of degree: 
(p-l)[N] + N-l + e N = 



pH - 1 for N odd 
p(N - 1) + 1 for N even 



in A. So for N odd to completely characterize it we have to fix its value in pH distinguished points and 
we are free to chose these points coinciding with the zeros of the operator B(A). For N even, we have to 
add to the B-zeros the values at the infinity, so that by using the corresponding interpolation formula for 

B(AgP- 1 )B(AgP~ 2 ) • • • B(Ag)A(A), we derive: 

B- 1 (A)A(A) = ^( — + ^:V4EE a ^ r» • ^ 

From the previous formula and the representation (I6.24t - (I6.25I >. we have that the local operators u^ and 
Pk,n are linear combinations of the monomials [J n rj N eNh (^^j ®ai,hi " " " Oa^^U" 1 for s < p, a\ € 
{1, [N]} and h, hi € {0, ...,p — 1}. For any monomial O ai ,h x • • • Ca s ,h s we can use the commutation 
rules (16.181 ) to rewrite it in a way that operators with the same index a are adjacent, we can order them in a 
way that a, < aj for i < j € {1, [N]} and we can apply the rule (16.151 ) to say if the monomial is zero or 
not. Finally, by using the property (I6.16I ). we have: 

•A(Zg) ( a ) 

llb^a,b=iy Z a/ Z b ~ ^b/^a) 

and so it is clear that all the non-zero monomials O aij hj • • • Oa s ,h 3 can be written in the form (16-2 1 b - □ 



6.2.2 Form factors of elementary operators 

As anticipated the interest in the above definition of elementary operators is the simplicity of their form 
factors: 

Lemma 6.3. Let {tk\ and \t' k ,) be two eigenstates of the transfer matrix T(A), then it holds: 

je N e N h k' 

/+, lc(ai.-.ar) i./ \ _ k,k>+h1 f , , /|i n (e N /io,{a},{a})m 

W C (hM)^Ml,hl),---,(ar,hrV Z k') - (0)he Nrr N c e N h Wo,M,M) m <™_ \ 1 1 U a,fe 

(6.29) 

where \\0^ h °>{ a }>{ a "\\ is the ([N] + rp — g) x ([N] + rp — g) matrix of elements: 

s (i" 1 ^) 4 ^ /ori ™ e {0 '-' p - me ^--^ (6 - 30) 

oStSSw ^ Cl(^o +9 )/ 2 /- ^ [N] - »•}, ff = E°*' ( 6 - 31 > 
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for any a € {1, [N] + rp — g}. Here, we have defined {b\, ...,b\N]- r } = {!> •••> [N]}\{oi, a r } with 
elements ordered by bi < bj for i < j and 

( (hi)\ e N''0+<*i([ N ]-'") 

f(. N ho,{a},{o}) = " ; (h . +tt ._ h) 7^ 1^ l^T 

Tir r-raj-1 tti-1 /^aj ' a j \ TT r f /a J 1 

llj=l ll/i=0 llj=U (h,0 (hi+^-h) Hlj=i+ll (h j+ h) (h,) J 

(-l)ELi(°i-i) H[ =1 g -([N]-r)a < (a < -l)/2 y((77 gi) ) 2 ) ^ (^0)2) 

x ni=i nJ N Jr - ^(faS ) 2 , (^ i+p ~ ai) ) 2 , .... rf^ 2 ..., (^ +p ^ } ) 2 )) ' 

(6.32) 

where V(xi, ...,£n) = r[i<&<a<N( x a ~~ x b) iS ^ g Vandermonde determinant. 

i e N^0 

'/A 

-e N /i | 

„h (k-h) \ SN P fjk-h)h N \ eN 



Proof. The operator ^ N Sn/i ( |j ) act in the following way on the state (t^ \ : 

(tk\VN 



(J)«n N thai ^ \ p 1 ' 2 i 

/ ( fc N)| 

x n (tf- ) ) 2 -tf >) ) 2 ) ™ N1 '-'" ( i ( ) ■ <6 - 33) 

!<6<a<[N] 116=1^61% J 



From the formula (16.191 ), it follows: 



n [N] 

11^,6=1 



Wl ,-, ! -,f/o i .->'/ N l^oi.hi n [N] TTai-l.^-^/Cfeb) (fci)/ (^-ftjx • 



So we can compute also the action of O ^ ■ ■ ■ O a ^ just taking into account the order of the operators in 
the monomial which leads by the scalar product formula to: 



(**l f ( ( S)"?(ai,fci),..,(a r ,h r )l*fc') - ( / (0K hTT N thoj E K p n^ a) ) h0eN 

\l% J 1 lj=l Sj / fei,...,fc N =l \ y / a=l 



X 



n 



? h (fe-/i) \ SN ^ /„[(A:-/i)-fc']fcN\ eN 

, °)^/irr N / 

r a<-l ttt / (hi-fe) / (ftj) _ (/y) / (Ai-h)-v-l 
i r t r 1 lj=i+i vhi I'laj >laj / 'hi ) 

1111 (Jhi-h) , (hj-aj) (hj-aj) , (/n-ft)\ 

i=l h=0 llj=iwai /^i ~~ ^aj /'/ai J 

11 , (^). 11 / (hi)N 

j=l W 6j (%. ) »=1 ^Woi J 

(W- (6-35) 
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Let us remark that the sum Y%i fe N =i reduces to 5^ N k , +h times the sum Y^k b k b =1 f° r trie presence 
of the ni=i 5k a Now we multiply each term of the sum by: 



i=nnii ri ^t^?-^?) 



e=±l i=l j=l h=—p+OLi 



(6.36) 



here the power +1 leads to the construction of the Vandermonde determinant: 

v((v^?, (e +p ~ ai) ) 2 , (i M ) 2 , (e +p - ar) ) 2 > (% i} ) 2 , .... (Cr; r) ) 2 )' ( 6 - 37 ) 



p-ai+l columns p-a r +l columns [N] -r columns 

and the sum V? , becomes sum over columns and after some algebra we get our formula. □ 

It is interesting to point out that the last [N] — r columns of the matrix ||o^ /l0 '^ a ^ a ^|| are just those of 
the scalar product. 

Remark 7. For the similarity of the model and representations considered, it is natural to cite the series 

of works lll43Ml46ll . There, in the framework of cyclic SOV-representations, first results on the matrix 
elements of local operators appear. However, it is worth saying that these quantities are there computed only 
for the restrictiorO of the T2-model to the generalized Ising model. In particular, the matrix elements of ui 
are computed and the results are not presented in a determinant form. 



7 Conclusion and outlook 
7.1 Results 

In this article we have considered the lattice sine-Gordon model in cyclic representations and we have solved 
in this case two fundamental problems for the computation of matrix elements of local operators: 

• Scalar products: determinant of N x N matrices whose matrix elements are sums over the spectrum 
of each quantum separate variable of the product of the coefficients of states, this being for all the 
left/right separate states in the SOV-basis. 

• Inverse problem solution: reconstruction of all local operators in terms of standard Sklyanin's quantum 
separate variables. 

Further, we have shown how these results lead to the computation of matrix elements of all local operators. 
At first, standard^! Sklyanin's quantum separate variables are suitable for solving the transfer matrix spectral 

17 It can be compared to the restriction to the case q 2 = 1 for the even sine-Gordon chain. 
18 Coinciding with the operator-zeros of one of the Yang-Baxter algebra generators, like B(X) or C(A). 
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problem. Indeed, the transfer matrix spectrum (eigenvalues & eigenstates) admits a simple and complete 
characterization in terms of Baxter-equation solutions in this SOV-basis. Then the inverse problem solution 
allows to write the action of any local operator on transfer matrix eigenstates as finite sums of separate states 
in the SOV-basis. Hence, the matrix elements of any local operator are written as finite sums of determinants 
of the resulting scalar product formulae. 

We have explicitly developed this program characterizing the matrix elements of the local operators u„ and 
a n by one determinant formulae in terms of matrices obtained by modifying a single row in the scalar 
product matrices. Moreover, we have constructed an operator basis whose matrix elements are in turn 
written by one determinant formulae. The matrices involved have rows which coincide with those of the 
scalar product matrix or with those of the Vandermonde matrix computed in the spectrum of the separate 
variables. 



7.2 Comparison with previous SOV-results 

In the literature of quantum integrable models there exist several results on matrix elements of local operators 
which can be traced back to applications of separation of variable methods. In this section, we try to recall 
the most relevant ones as they allow for an explicit comparison with our results. It leads to a universal 
picture emerging in the characterization of matrix elements by SOV-methods. 



7.2.1 On the reconstruction of local operators 

One important motivation for our work was to introduce a well defined setup which allows to solve the 
longstanding problem of the identification of local operators in the continuum sine-Gordon model thanks to 
the reconstructions achieved on the lattice. Then, it should also allow for the identification of form factor 
solutions of the continuous theory by implementing well defined limits from our lattice formulae. 

Even if methodologically different, it is worth recalling the semi-classical reconstruction presented by Ba- 
belon, Bernard and Smirnov in [93] for chiral local operators of the restricted sine-Gordon model (in the 
infinite volume) at the reflectionless points, /3 2 = 1/(1 + v) with v E Z- . The classical sine-Gordon model 
admits a SOV description: each ?7,-soliton solution ip(x, t) of the equation of motion can be represented in 
terms of n-separate variables Aj, which in the BBS choice Ir93ll lead to the representations: 

n . 

3=1 3 



where the Bj are integrals of motion. The formula (T7.U represents classically a SOV reconstruction of the 
local fields when restricted to the n-soliton sector. In (93|], this reconstruction has been extended to the 



quantum model in each n-solitons sector by quantizing the separate variables^] A, and the conjugate mo- 



menta as operators which generate n independent Weyl algebras with parameter q = e^ 1 -? 2 = e l «. This 
extension and the consequent identifications of primary fields and their chiral descendants in the perturbed 



>9 A fundamental point is the introduction of appropriate Hermitian conjugation properties and the characterization of the spec- 
trum of the Weyl algebra generators. 
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minimal models Mi.i +P , are justified by the following indirect but strong arguments: a) The ra-multiple 
integrals of the form (36) in H93I1 which represent the n-solitons to n-solitons form factor^ of chiral left 
operators at the reflectionless points are reproduced from the semi-classical limit, b) The counting of these 



form factor solutions allows the reconstruction of the chiral characters of Mi i+j, 19411 . Further support to 



7.11) was given by Smirnov's work on semi-classical form factors 21 ! of the continuous KdV model in finite 



volume 119711 : there the form factors of Il93ll were reproduced by taking the infinite volume limit of the KdV 



semi-classical ones. 

Let us remark that in our lattice regularization of the sine-Gordon model, choosing as quantum separate 
variables standard 22 ! Sklyanin's ones, the reconstruction of the exponential fields is not of the simple form 
given in (17.11) . Then, the following question is relevant: is it possible to find a SOV representation of 
the quantum lattice sine-Gordon model where the exponential fields are simply written as a product of 
generators of the SOV representations? 

A natural idea can be to implement a change of basis in the quantum separate variables from Sklyanin's 



ones to a new set; an interesting example of this approach was used by Babelon in M99H which has provided 
a simple reconstruction of the lattice quantum Toda local operators in terms of a set of quantum separate 
variables defined by a change of variables from the Sklyanin's ones. However, it is worth pointing out that, 
for a new SOV representation to be really useful for the computation of matrix elements, it should not only 
give a simple reconstruction of the local operators but also keep the solution of the transfer matrix spectral 
problem and the scalar product formulae as simple as for original Sklyanin's variables. Let us comment that 
a reconstruction like (17.11 ) can be formally derived at the quantum level implementing the special limit 

Kn/i -> +00 (7.2) 
on the following reconstruction formula of the lattice sine-Gordon model: 

{q- 1 ^ + k\ ) 



l\/2 t-- 2 
v n rv n 

The result for an even chain reads 



(q "11" a 



U n A- 1 (/in,-)B(^ n ,_)B- 1 (^ ni+ )A( Aini+ )U- 1 . (7.3) 



-2ft 



e 



2h N-1 



n N II ^U- 1 , h£{l,.., P -l} (7.4) 

1 la^n,a=l ?a a=l 

in terms of Sklyanin's separate variables. It is possible to argue that the previous limit can be consist ently 
interpreted as a chiral deformation of the lattice sine-Gordon model to chiral KdV models HQ- In 
a future publication, we will analyze the cyclic representation^] of these chiral KdV models showing our 
statement on the reconstruction formulae and computing the matrix elements. It is worth pointing out that 
the form factors of the r.h.s. of (I7.4t are trivial to compute in our SOV framework and are written by one 
determinant formulae which differ w.r.t. the scalar product only for/iG {!,.., p — 1} rows. 



20 Note that these are solutions of the form factor equations and so they surely represent local fields in the S-matrix formulation 
of the restricted sine-Gordon model. 

21 On the basis of this last Smirnov's work, Lukyanov has introduced his conjecture for the finite temperature expectations values 
of exponential fields in finite volume for the shG-model [98]. 

22 The same is true if we take the products of the operator zeros of C(A) but also of A(A) and D(A), i.e. for all possibilities to 
construct the SOV representations by the simplest Sklyanin's method. 

23 Here, we are referring to compact representations of chiral KdV models where the generators of the local Weyl algebras are 



unitary operators. The spectrum of the non-compact versions was instead analyzed by SOV and Q-operator method in [ 102]. 
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7.2.2 On the matrix elements of local operators 



In the case of the quantum integrable Toda chain 18811 . Smirnov 19711 has derived in the framework of 



Sklyanin's SOV determinant formulae for the matrix elements of a conjectured basis of local operators 
which look very similar to our formulae. The main difference is due to the different nature of the spectrum 
of the quantum separate variables in the two models. In fact, in the case of the lattice Toda model, Sklyanin's 
measure is continuous (continuous SOV-spectrum) while it is discrete in the case of the cyclic lattice sine- 
Gordon model. The elements of the matrices whose determinants give the form factor formulae are then 
expressed as "convolutions", over the spectrum of the separate variables, of Baxter equations solutions plus 
contributions coming from the local operators. In the case of Smirnov's formulae they are true integrals, the 
SOV-spectrum being continuous, while in our formulae they are "discrete convolutions" , the SOV-spectrum 



being discrete. Let us comment that the need to conjecturd I the form of a basis of local operators in 119711 is 
due to the lack of a direct reconstruction of local operators in terms of Sklyanin's separate variables. 



In the case of the infinite volume quantum sine-Gordon field theory, the form factors of local operators |67 1 
have also a form similar to the one predicted by SOV. This similarity can be made explicit considering the 
n-soliton form factors for the restricted sine-Gordon theory at the reflectionless points in formula (31) of 



09311 . Then, for the local fields interpreted as primary operators in [93], the corresponding form factors can 
be easily rewritten as determinants ofuxn matrices whose elements are integral convolutions of n-soliton 
wave functions (the ^-functions (32)) plus contributions coming from the local operators. 

The rough picture that seems to emerge is that by performing the IR limit on our lattice form factors the 
lattice wave functions factorized in terms of Q-operator eigenvalues have to converge to the infinite volume 
n-soliton wave functions. To which extent this picture can be confirmed and clarified by a detailed analysis 
of the thermodynamic limit starting from our lattice sine-Gordon model results is of course an interesting 
question to which we would like to answer in the future. 



7.3 Outlook 

It is worth mentioning that we didn't succeed yet to express the matrix elements of discretized exponential 
of the sine-Gordon field in terms of one simple determinant formula. Hence our next natural project is 
the simplification of the present representation; this is also important in view of the attempt to extend 
our results from the lattice to the continuous finite and infinite volume limits. The main goal here is to 
derive the known form factors of the IR limit, starting from our lattice form factors, in this way solving 
the longstanding problem of the identifications of local fields in the S -matrix characterization of the infinite 
volume sine-Gordon model. 

Beyond the sine-Gordon model we want to point out the potential generality of the method we have in- 
troduced here to compute matrix elements of local operators for quantum integrable models. The main 
ingredients used to develop it are the reconstruction of local operators in Sklyanin's SOV representations 
and the scalar product formulae for the transfer matrix eigenstates (and general separate states). The emerg- 
ing picture is the possibility to apply this method to a whole class of integrable quantum model which were 

24 The consistency of this conjecture is there verified by a counting argument based on the existence of an appropriate set of null 
conditions for the "integral convolutions". 
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not tractable with other methods. This is in particular the case for lattice integrable quantum models to 
which the algebraic Bethe ansatz does not apply. The first remarkable case is given by the T2-model in 
general representations which are of special interest for their connection to the chiral Potts model. This will 
be the next model that we will analyze by our technique due to the similarity of its cyclic representations 
with those of the sine-Gordon model. 

There are also many other examples which are interesting and for which, on the one hand, the reconstruction 



of the local operators can be deduced from 02911 and, on the other hand, the description of the spectrum can 
be given by Sklyanin's quantum separation of variables. For all these models the possibility to apply our 
method for the computation of matrix elements is very concrete and moreover the results are expected to 
have a completely similar form to the ones shown in the present article. 
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